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1 Introduction 



Remarkable relation between the spectral theory of Schrodinger equation [1] 
and vacuum eigenvalues of the so-called Q operators of Conformal Field 
Theory (CFT) was discovered a few years ago in [2]. The Q-operators were 
introduced in [3] as a mathematical tool to describe integrable (Yang-Baxter) 
structure of c < 1 CFT, where they play the role of the Baxter's Q-matrix 
[4]. It was observed in [2] that for several degenerate Virasoro module the 
vacuum eigenvalues of the Q-operators essentially coincide with the spectral 

determinants of the Schrodinger operator — -^-^ + V(x) with the potential 
V{x) = x 2a , where the exponent a is determined by the Virasoro central 
charge c. This relation was extended in [5] to all Virasoro modules with 
c < 1. On the other hand, the Virasoro vacuum is one among a sequence of 
eigenstates of the operators Q in the whole Virasoro module. The question 
arises if such relation to the Schrodinger operator exists for the higher-level 
eigenvalues. 

In this note we answer this question positively and identify the Schrodinger 
operators associated with the higher-level eigenvalues. For an eigenstate on 
the level L of the highest-weight Virasoro module the potential has the form 

V{X) = fcll + J* _ 2 |L £ log(x2 «+2 _ Zk) ^ (1) 

fc=l 

where a and £ are related to the central charge c and the highest weight A 
of the Virasoro module as follows 

c-1-isL, A = ^+i> 2 -^ 2 . (2) 

a + 1 16 (a + 1) w 

The Zk, k = 1,2, . . . L are pairwise different nonzero complex numbers sat- 
isfying the system of L algebraic equations 



3 = 1 



\ - Zk{4 + (3 + a)(l + 2a)z k Zj + a(l + 2a)z]) aZk _ n 



with k = 1, 2, ... L. In what follows we refer to the potentials (1) as the "Q- 
potentials" . The arguments presented below strongly suggest that appropri- 

j2 

ately defined spectral determinants of the Schrodinger operator — -P-j + V r (a;) 
with the Q-potentials V(x) coincide, up to certain trivial factors, with the 
eigenvalues of the Q-operators. 



The Schrodinger equation with the monstrous potentials (1) are unlikely 
to be of much direct interest in quantum mechanics. On the other hand, 
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the higher-level eigenvalues of the Q-operators are definitely of interest in 
integrable quantum field theories. One of the reasons is their relation [3,6] to 
the boundary states associated with certain integrable boundary flows [7,8]. 
This work was motivated by possible applications to this kind of problems. 

We did not attempt at making this paper self-contained. The definitions 
of the Q-operators and the local and nonlocal integrals of motion are as in 
Ref. [3]. The arguments in Section 3 are meant to be read in conjunction 
with Refs. [2,5]. 



2 The Q-operators 

Here we only list some basic properties of the Q-operators used in the dis- 
cussion below. We refer the reader to [3, 9] for explicit construction of the 
Q-operators and more details. 

The Q-operators are actually the operator- valued functions Q±(s) of a 
complex parameter s 1 . The operators Q±(s) act in the Virasoro module 
Va- It will be convenient to parametrize its highest weight A and the Vi- 
rasoro central charge c < 1 as in Eqs.(2), where we assume that the real 
parameters I and a are in the range i > —1/2, a > 0. The operators Q±(s) 
form a commutative family and they also commute with the zero-mode Vi- 
rasoro generator Lo- The Virasoro module Va admits the standard level 
decomposition, 



Va = ©?=o vi L) , L vi L) = (A + L) V { £ ] , (4) 

3vel 

Q-operators, 



and the level subspaces V/^ are invariant with respect to the action of the 



Q± : V1 L) - vi L) . (5) 
In particular, the Virasoro vacuum | A ) is an eigenvector of the Q-operators, 

Q±(s) \A) = Q ( r c) ( S ) | A). (6) 

It is the vacuum eigenvalues Q± ac \s) that have received most of attention 
in the Refs. [2,5]. Here we will discuss generic eigenvalues Q±(s) of the 

lr The operators Q±(s) discussed here coincide with Q±(A) of [3,9] with 

3 (lfJ S - 

The parameter a is related to (3 used in [3,9] as 1 = 



A 2 = 



r- 2 ' 
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operators Q±(s) in Va- As the functions of s they enjoy the following basic 
properties [3,9] 2 ; 

2£-\-l 

(i) The functions A±(s) = (— s) T ^~ Q±(s) are entire functions of s, and 

A±(0) = 1 . (7) 

(ii) The eigenvalues Q±(s) have the leading asymptotic behavior: 

logQ ± ( S ) = — * (_ a )^+0(l) (8) 

at |s| — > cx3 , arg(— s) < n. 

(Hi) The functions A±(s) satisfy the functional relation (the "quantum 
Wronskian condition") 

q^A + {qs)A^{q- 1 s) - q~^ A^q' 1 s)A_(qs) = q £+ ^ - , (9) 

where 

q = e Th. (10) 

The equation (9) follows directly from the construction of the Q-operators 
[3,9]. The statements of analyticity, (i) and (ii) above, and especially the 
domain of validity of the asymptotic form (8) , presently have somewhat less 
solid status. For the vacuum eigenvalues Q± ac \s) these properties follow 
from the results of [5]. The solution at c = —2, where all eigenvalues can be 
found explicitly [3], also supports the above analyticity statement (i), (ii) 
(see footnote 2 on this page) . But it is fair to say that a solid proof in a general 
case is still missing. Nonetheless, in the further discussion we assume these 
properties. 

Entire functions A±(s) having the asymptotic behavior (8) are charac- 
terized by positions of their respective zeros s^ in the complex s-plane. The 
zeros accumulate towards s = oo along the positive real axis. The quantum 
Wronskian equation (9) then leads to infinite set of the Bethe Ansatz equa- 
tions for the numbers s^. The analysis of these equations in Appendix A 
shows that there is a discrete (albeit infinite) set of functions Q±(s) that 
satisfy all the conditions (i — Hi) above. Clearly, all eigenvalues of the Q- 
operators can be found among this set. But there is much evidence that the 

2 In writing Eqs.(8),(9) below we assume that (1 + a)/ (2a) is not an integer. If it equals 
to N = 1, 2, . . ., the asymptotic form (8) changes to s N log(— s) and (9) must be modified 
accordingly. The free-fermion point a = 1 (i.e. c = —2) is one of these "resonance" cases. 
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converse is also true — any pair of functions Q±(s) satisfying {i — Hi) are 
eigenvalues of the pair of operators Q±(s) in Va- Given a solution of (i — iii), 
the level L of the associated eigenvector can be identified as follows. Correc- 
tions to the leading asymptotics (8) are expressed through the eigenvalues of 
the local and nonlocal integrals of motion, according to the expansion (26) 
below. Since the first local integral Ii coincides with the Virasoro generator 
Lo — c/24, the first subleading terms to the asymptotics (8) have the form 
(a > 1/2): 

logQ±(s) ^ -4ttN^ + C 

where the constant term does not depend on s (but, of course, depends on 
a state from Va) and L is the level. The results of Appendix A suggest that 
for generic c and A the number of those solutions of (i — in) which also have 
the subleading behavior (11) with a given integer L > equals the number 
of integer partitions of L. It thus coincides with the dimension of the level 
subspace V^ of a generic non-degenerate Virasoro module Va- 



3 Schrodinger equation 

It is known (see e.g. [1,5]) that appropriately defined spectral determinants 
D±(E) (Eq.(21) below) for the Schrodinger equation 



-— T *(x) + V(x) ^(x) = E^/(x) 
dar 



with the potential 



V ( x ) = V (vac) {x) = W+}) +x 2a 



(12) 



(13) 



obey all the above properties (i — Hi), if one makes the identifications 

D±{v 8 ) = A ± (s) , (14) 

where 

2 v^r(| + ^; 



2a 
1+a 



(15) 



Thus D±(E) are related to some eigenvalues of the Q-operators. In fact, 
the spectral determinants of (12) with V(x) = V^ vac \x) correspond to the 



716 



Higher-level eigenvalues of Q- operators and Schrddinger equation 



vacuum eigenvalues Q± . This correspondence in the case £ = is the 
original observation of [2]. It was proven in [5] by combining Eq.(9) with 
known analyticity of Q± ac ^ as functions of the parameter I. 

It is important that the properties (i — Hi) of the spectral determinants 
can be derived using only the following three properties of the potential (13) 
as the function of complex x: 

• (I) Symmetry: 

V(qx)=q- 2 V(x) , (16) 
where q is defined in (10). 

• (II) Asymptotic behavior: 

V(x) -> e ^^ 2 ^ + o(l) as x^O, (17) 
V(x) -► x 2a + o(l) as x -> oo . (18) 

• (III) Monodromy properties: For any value of E all solutions of (12) 
with this potential are monodromy free everywhere except for x = 
and x = oo. 



Here we observe that V i - vac \x) is not the only potential which meets all of 
the conditions (I — III). We show in Appendix B that the most general 
potential V(x) with these properties is given by (1), with the parameters 
Zk satisfying the Eqs.(3). Note that for L > the potential (1) has poles 
away from zero, at x 2a+2 = zi, Z2, ■ ■ ■ , zl- The equations (3) then guarantee 
that these poles do not violate (III) ■ Thus (I — III) are the characteristic 
properties of all Q-potentials. 

Given a Q-potential, Eqs.(l),(3), the associated spectral determinants 
D±(E) are defined as follows. For %tel > (1 + a) L - 3/2 the equation (12) 
has a solution i/j(x,£,E) uniquely specified by the condition 

ip(x,EJ): ip(x,E,£) = x i+1 + 0(x £+3 ) as x^O. (19) 

This solution can be analytically continued outside the domain ^ke£ > (1 + 
a) L — 3/2. Let tp(x,E, —£ — 1) be another solution of the same equation 
(12) obtained from ip(x,E,£) by appropriate analytic continuation in £ 3 . 

3 In general, analytic continuation of ip(x,E,£) in I yields a multivalued function of t, 
and tjj(x, E, — i— 1) at different branches solves (12) with different potentials, corresponding 
to different solutions of the same equations (3). However, since I enters the Eq.(3) only 
through A, Eq.(2) , there exists a continuation ^ — > — ^ — 1 which leaves the potential (1) 
unchanged. 



V.V. Bazhanov, S.L. Lukyanov and A.B. Zamolodchikov 



717 



For generic values of I the solutions ip+(x,E) = ip(x,E,£) and if)-(x,E) = 
ip + (x,E, —t — 1) are linearly independent, since their Wronskian 

W[ip + ,ip_] = i; + d x ip_ - y>- d x ip + 

equals to -{21 + 1). 

For certain isolated values of E one of these solutions decays at x — > 
+oo. Let {E^}°^ =1 and {E^}^ =1 be ordered spectral sets defined by the 
conditions 

V+(x)^V + (z,£ fc + )^0, Vfe (x) = ^-(x,E-)^0, (20) 
as x — > +oo. One then defines 

^(^)=n( 1 -^) • ( 2i ) 

k=l k 

Straightforward WKB analysis of the equation (12) shows that 

E k= v {k-\±^)^ (l + o(l)) , as fc^oc, (22) 

and therefore at a > 1 these products converge for all finite E. For < a < 1 
the definition (21) has to be modified by including the Weierstrass factors, 
but this does not affect our arguments. 

The spectral determinant D±{E) thus defined are entire functions of E. 
Asymptotic behavior of these functions at large E follows from (22), 

logD ± (E) = (-E/v) 1 -^T^log(-E) + 0(l), arg(-£)<7r. (23) 

And one can repeat the arguments of [5] to show that the spectral deter- 
minants (21) satisfy exactly the same functional relation (9), as the direct 
consequence of the properties (I — III) of the Q-potential. We conclude that 
the pair of functions D±(E) satisfies, upon the identifications (14), all the 
analytic conditions (i — Hi) of Section 2. This shows that every Q-potential 

(I) corresponds to some eigenvector of the Q-operators, and the associated 
eigenvalues Q±(s) coincide, up to the factors (— s)^ i , with the spectral 
determinants D±(vs). Moreover, using the WKB expansion of the solutions 
of (12), it is possible to evaluate corrections to the leading asymptotics (23) 
for any Q-potential (1). For the first such correction one obtains exactly 

(II) with L standing for the number of nonzero poles z±, Z2 , . . . ,zl of (1). 
This confirms the identification of L in (1) with the level. 

At this point the question arises: is this description of the spectrum 
of the Q-operators complete? In other words, is the family of potentials 
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(1) large enough to accommodate all the eigenvalues Q±(s)? Answering 
this question requires enumerating all solutions of the algebraic equations 
(3) for an arbitrary level L. We did not accomplish this task. However, 
an elementary analysis of a few lowest levels suggests positive answer to 
this question. For instance, for A / the level 1 subspace is one- 
dimensional. Correspondingly, for L = 1 the equation (3) has only one 
solution 



Zl = Z = ^ 1 + a) A. (24) 
a 



Therefore the spectral determinants D±(E) of (12) with 



V{x) = x- 2 | x 2a+2 + £(£ + 1) + 4a + 4 



4(a + l)(2a + 3)z 8(a + 1) 2 z 2 , 

x 2a+2 _ z ( x 2a+2 _ -,)2 ( 



coincide with the eigenvalues of Q±(s) corresponding to the eigenvector 
L_i | A). Note that z turns to zero at A = 0, in agreement with the fact 
that in this case L_i | A ) becomes a null-vector. 

For L = 2 the system (3) is solved by elementary methods. For generic c 
and A it has exactly two solutions (of course we do not distinguish solutions 
which differ by permutations of Zk). We can recall that for a non-degenerate 

(2) 

Virasoro module the dimension of equals 2. Again, one of the solutions 
disappears (in the sense that one or both of the numbers zi, Z2 turn to zero) 
in three cases A = 0, A = and A = as expected, because at 

these values of A null- vectors appear at the level 1 or 2, and the dimension 

(2) 

of drops down. Similar situation is observed at the level 3. 



4 Local integrals of motion 



Finally, let us make remark on the eigenvalues of the local Integrals of Mo- 
tion (IM) and their relation to the semiclassical analysis of the Schrodinger 
equation with the Q-potentials (1). Asymptotic expansions of the opera- 
tors Q±(s) around the essential singularity at s = oo (in the sector — ir < 
arg(— s) < 7r) are expressed in terms of the infinite sequence of the local IM 
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I 2n -i and "dual nonlocal IM" H± [3], 



oo 

^ / n I* I xl+H ^ ST^ l , (l-2n)(l+g) 

IogQ±(a) ~ j— (-a) a- +C + 2^(-i/s) 2- M2»-i 

COS V 2 ai n= i 

00 



n=l 



where the coefficients b n and c ra have explicit form 



b - ( - 4) (i 1 arrf 1 - 2re vr (2n - 1)(1 + a) i 

c n = (-I)™" 1 4 na {l + a)- 2n . (27) 



The asymptotic expansions of the eigenvalues Q±(s) are thus expressed 
through the eigenvalues of the local and nonlocal IM. Note that (26) has 
exactly the structure which appears when one performs the WKB expansion 
of the solutions of the Schrodinger equation (12) with a Q-potential (1). The 
correspondence (14) suggested in the previous section implies that the coef- 
ficients generated by this WKB expansion must agree with the eigenvalues 
of the operators l2n-i an d H^. 

Whereas the eigenvalues of the dual nonlocal IM are relatively difficult 
(but not impossible) to calculate, the local IM \2n-1 have explicit expres- 
sions in terms of the Virasoro generators L n (see Eq.(ll) of Ref. [6]), and 
their eigenvalues in any given level subspace can be found by a di- 
rect calculation. Thus, for the vacuum eigenvalues of the first few l2n-i we 
have [6] 



& ac \A) = A- — , 
1 v ) 24 , 

, r > (A) _ A2 _£±? A + £(|i_M, (28) 

Avac) (A .3 c + 4 2 (c + 2)(3c + 20) 

h (Aj " A 8 + 576 
c(3c+ 14)(7c + 68) 
290304 ' 



where the Virasoro central charge c is related to the parameter a as in Eq.(2). 
We have checked by explicit diagonalization at the levels L = 1, 2, 3 that the 
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eigenvalues of the first three local IM agree with the following expressions 
h = l{ vac \A + L), 

h = 4 mc) (A + L) + ^— z k > (29) 

k=i 

t r(™ c )^A r x , / A + L — 1 3 + Q \ 

h = * (A+L) + (^T^ + 24(rT^) x 

3a (3 + 2a) s~\ 3 a 2 2 

l + « ^ Zfe + 2(5 + 2a)(l + a) 2 ^ Zfc ' 

in terms of the same parameters Zk that enter (1). One can verify that the 
same expressions come out from the WKB analysis of Schrodinger equation 
(12),(1). 
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Appendix A 

In this Appendix we assume that a and £ take some generic values so that the 

Virasoro module Va is non-degenerate. Moreover we assume that a > 1 4 . 

2^+1 

Then the entire functions A±(s) = (— s) T 4 Q±(s) are uniquely determined 
by their leading asymptotics at s — > 00 (8), the normalization (7), and the 
locations of their respective zeroes Sk, {k = 1, 2, . . .) in the complex s-plane. 

4 For < a < 1 the definition (A. 3), (A. 4) below may need modification, but the state- 
ment about the one to one correspondence between the eigenstates and the sets (A. 6), (A. 7) 
most likely remains valid. 
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From the quantum Wronskian relation (9) it follows that for any eigen- 
vector of Q±(s) the zeroes of A±(s) satisfy the same infinite system of 
the Bethe Ansatz equations 



A±(s 



± „2 



g^+D . (A.l) 



The eigenvalues corresponding to different eigenvectors are distinguished by 
a specific phase assignment in this equation 



2^i l0g 



A±(s±q 2 )~\ 21 + 1 \ 



^WTa)- mt+ 2> (A - 2) 



determined by a choice of integers {w-^}^ =1 . These integers, of course, 
depend on the choice of branches of the logarithm in the left hand side of 
(A. 2). However, once these branches are appropriately fixed (see below), 
every eigenvalue is characterized by a unique set of the integers {m^}. 

As follows from (ii) in the Section 2, the zeroes accumulate toward 
infinity along the positive real axis. In fact, explicit solutions for special 
values of c, as well as numerical analysis, give all support to the following 



Conjecture: For all eigenvalues at the level L, and for real t > (1 + 
a)L — 3/2 all zeroes of A + (s) are simple, real and positive. The zeroes 
of A_(s) enjoy the same property for I < — (1 + a)L + 3/2. 



Assuming this property one can unambiguously fix the phase of the 
logarithm in (A. 2). Consider, for example, the zeroes {s^} of A + (s) for 
£ > (1 + a)L - 3/2 and define the L.H.S. of Eq.(A.2) as 

1 °° 

L.H.S. of Eq.(A.2) = - V arg(l - q 2 s k / Sj ) , (A.3) 

where for real positive x 

-ir < arg(l - q 2 x) < +ir . (A.4) 

Since all zeroes are assumed to be distinct and positive they can be 
uniquely ordered by their value < < < . . . . Then for each 
eigenvalue A + (s) all the integers will be distinct and uniquely defined. 
For other values of £ the L.H.S. of (A. 2) is defined by means of the analytic 
continuation of (A.3) in I. 
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Obviously, not every set of integers {m k } corresponds to an eigenstate 
of Q±(s) in Va- Indeed, it follows from (8) that 

4 = (fc-i±^±±)^ (l + o(l)) , as k^oo. (A.5) 

Hence (A. 2) implies that for a given eigenvalue the sequence of integers 
stabilize at large k, i.e. = k for sufficiently large k. The sets {m^} associ- 
ated with different eigenvalues differ in finitely many first entries. Therefore 
the most general pattern in the set {m^} (or {m^}) can be obtained from 
the vacuum set (m^ = k, for all k = 1, 2, ... , oo) by deleting a certain num- 
ber of entries (we denote this number by M) and adding the same number 
of distinct non-positive integer entries. It can be written as 

m + _ i 1 - n M-k+i > for k = 1, . . . M , 
k ~ I N k _ M {n+) , for k > M + 1 . 1 ° J 

Here n 1 * 1 = {nf,nf,...,rij^} denote two increasing sequences of positive 
integers 1 < nf < n\ < . . . < ri^ and 1 < < n% < ■ ■ ■ < n^, 
with M > 0; and Nj(n + ), j = 1,2, ...,oo, denotes j-th. element of the 
increasing sequence of consecutive positive integers with deleted entries nf, 
i = 1, . . . M: 

{tf(n+)} = { 1,2,... ,#,...}. (A.7) 

Using the standard technique based on the Destri-de Vega equation [10,11] 
associated with (A. 2) and (A. 6) one can easily calculate the next subleading 
term in the large s asymptotics of the corresponding eigenvalues A + (s). The 
result is given by (11) with 

M 

^ = EK ++n r- 1 ) • ( A - 8 ) 

For a given value of L there are exactly p(L) sets sjf^ = {n + \n~}, with 
k = 1,2, ...p(L), satisfying (A.8) 5 , where p(L) stands for the number of 
integer partitions of L. These sets naturally enumerate all eigenvectors in 
the level subspace 6 . As for the set {m^} describing the zeroes of 
A-(s), it is obtained from {m^} by interchanging the roles of n + and n~ in 
Eq.(A.6). 

5 The counting follows from the following elementary partition identity 

oo oo n oo n 2 n 

E p(») * n = E Tir = E - (*)» = II (! - ^) ( i-i < !) • 

n = n=0 yZ)n n=0 Wn fc=l 



For example, at L = 4, there are five sets <S£ \ {4|1}, {1|4}, {2|3}, {3|2}, {1,2|1,2}. 
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As explained in [5] the formula (26) allows one to derive the asymptotic 
expansions for the zeroes of the eigenvalues A±(s). Indeed, let hn-i an d H„ 
are the eigenvalues of the local IM and the dual nonlocal IM corresponding 
to a certain eigenvector of Q±(s) determined by the set of phases (A. 6). 
Then substituting (26) into (A. 2) one obtains 

l+g JS, (l-2n)(l+o) J^2- , . ~ 

V* + E ("**) ~ a & 7 2n-l + E (^)~™ ( } 7n ^ 

n=l n=l 

^ mf-Ji^l, (A.9) 

where k = 1, 2, . . . and 

& = (-If vr- 1 cos(7r(n-l/2)/a) 6 n , 

7n = (-l) n+1 tt- 1 sin(Trna) c„ . (A.10) 

On the other hand the Bethe Ansatz equations (A. 2) complemented by the 
asymptotic condition (A. 5) can be solved numerically [12]. We have per- 
formed these calculations for the zeroes of A + (s) in a particular case a = 4 
with various choices of the parameter I and the integer phases {m^} as in 
(A. 6) above. For each of these choices we then calculated a few first coeffi- 
cients in the asymptotic expansion (A.9) and thus determined the eigenvalues 
hn-i numerically (Note, that for a = 4 the second sum in (A.9) vanishes.). 
The results are in very good agreement with those obtained from an explicit 
diagonalization of the first three local IM (29). The correspondence was 
verified for all eigenvectors at the level L < 5. 



Appendix B 



Here we explain why the form (1) describes the most general potential which 
satisfies the conditions (/ — 777) in Section 3. Indeed, let us write V(x) in 
the form 

V(x) = +x 2a + v(x) . (B.l) 

ar 

The conditions (/ — /i7) imply that 

v(x) =x~ 2 F(x 2a+2 ) , (B.2) 

where F(x 2a+2 ) is a rational function of its argument, vanishing at x = 
and bounded at x = oo, 

F(0) = 0, |F(oo)|<oo. (B.3) 
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Therefore, unless F = 0, the potential (B.l) has poles at finite nonzero values 
of x. Let xq 7^ be the position of any such pole. Consider the Laurent 
expansion of the potential (B.l) near this pole 

oo 

V{x) = (* " *o) k Vk ■ (B.4) 

k=r 

It is well known (see e.g. [13]) that coefficients of this expansion should be 
constrained as 

r>-2; y_ 2 = m(m + 1) ; V 2k -i = , k = 0, 1, . . . m , (B.5) 

where m = 0, 1, 2 . . ., to ensure that both linear independent solution of 
(12) are single valued functions x in the vicinity of xo for all values of E. 
One can check that for m > 2 there are too many conditions in (B.5) which 
cannot be simultaneously satisfied for generic values of a and I. So the only 
case to consider is m = 1 which requires 

y_ 2 = 2, V^ 1 = V 1 = 0, (B.6) 

for every pole of (B.l). These conditions, together with (B.2) and (B.3), 
lead to (1) and (3). 
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